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Full counting statistics (PCS) of electron transport is of fundamental importance for a deeper 
understanding of the underlying physical processes in quantum transport in nanoscale devices. Here 
we investigate the backaction of a charge detector in the form of a quantum point contact (QPC) 
on the FCS of a biased double quantum dot (DQD). We show that this inevitable QPC- induced 
backaction can have profound effects on the FCS under certain conditions, namely enhancing current 
flow through the DQD, and changing the shot noise from being sub-Poissonian to super-Poissonian. 
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I. INTRODUCTION 

Current fluctuations in nanoscale systems provide key 
insights into the nature of charge transfer beyond what is 
obtainable from a conductance measurement alone (see, 
e.g., Refs. [l] and[3 for recent reviews). An in-depth un- 
derstanding, however, may require us to go beyond the 
first-order and even the second-order current correlation 
functions (corresponding to the average current and shot 
noise respectively) to study the full counting statistics 
(FCS)^i^ which yields all zero- frequency correlation func- 
tions at once. Real-time detection of the tunneling of in- 
dividual electrons, an important step towards experimen- 
tal measurement of the FCS, has recently been achieved 
in various quantum-dot (QD) systemsi^"— In particular, 
since its measurement in a single QD for the first time,»S 
FCS has become an important experimental tool to ex- 
amine interaction and coherence effects in nanoscale sys- 
tems under out-of-equilibrium conditionsi^^— More re- 
cently, FCS was applied to characterize correlations in 
both classical and quantum systemsJ^ 

However, a pronounced effect known as backaction on 
the FCS of electron transport^ii^ is inevitably introduced 
during measurements made by even the most noninvasive 
detectors such as a quantum point contact (QPC) Very 
recently, such backaction has been investigated experi- 
mentally in a single QDJ^iii In contrast to a single QD, a 
double quantum dot (DQD)^ involves coherent coupling 
between two different dots, and therefore can be used 
to demonstrate prominent coherent effects. The FCS 
for DQDs has been studied theoretically^^ and experi- 
mentally only for noise properties In addition, for 
a zero-bias DQD, the effects of charge-detector-induced 
backaction was studied theoretically^ to explain exper- 
imental observations on inelastic electron tunnelingi^ 
However, to the best of our knowledge, the impacts of 
charge-detector-induced backaction on the FCS in these 
QD systems have not yet been studied. 



In this paper, we investigate the FCS of electron trans- 
port through a biased DQD under measurement by a 
charge detector. We demonstrate that this inevitable 
backaction can indeed have profound effects on the FCS 
under certain conditions. In particular, it enhances the 
current fiow through the DQD, and can change the na- 
ture of the shot noise from being sub-Poissonian to super- 
Poissonian. These QPC-backaction-induccd effects are 
expected to be experimentally observable with currently 
existing technologies. Apart from a deeper understand- 
ing of experimental observations, this study may also 
shed light on how to control quantum transport processes 
in these QD systems. Also, we show that when the mea- 
suring device is absent (i.e., when the charge detector is 
decoupled from the DQD), we recover the known results 
for both current^ and shot noisci^ 

This paper is organized as follows. In Sec.|TTl the model 
is explained and the formulas for calculating the FCS are 
derived. In Sec. lIIIl the backaction by the charge detector 
is analyzed. Summary is given in Sec. IIVI 



II. FULL COUNTING STATISTICS OF A 
DOUBLE QUANTUM DOT 

We focus on a setup consisting of a lateral DQD, 
which is coupled to the source and drain electrodes, and 
measured by a nearby QPC [see Fig. (Ua)]. The lat- 
eral DQD is formed by properly tuning the voltages ap- 
plied to the corresponding gates. Here we consider a 
Coulomb-blockade regime with strong intradot and in- 
terdot Coulomb interactions, so that only one electron is 
allowed in the DQD system. The states of the DQD are 
represented by the occupation states |1) and |2), denoting 
one electron in the left and the right dots, respectively 
[see Fig.mb)]. 

The total Hamiltonian of the whole system can be writ- 
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FIG. 1: (Color online) (a) Schematic diagram of a DQD cou- 
pled to two electrodes (S and D) via tunneling barriers. A 
QPC used for measuring the DQD electron states also yields 
backaction on the DQD. (b) Electronic transition between two 
eigenstates \g) and |e) of the DQD (with an energy difference 
A) can be induced by the charge detector QPC. The energy 
detuning e between the two single-dot levels (dashed lines) 
can be varied by tuning the gate voltages. 

ten as 

H = i?DQD + Hicads + Hqpc + Ht + -ffdot, (1) 

where (we set h= 1 ) 

s 

HdQD = -^CTz + ^CTx, (2) 

H\cads = '^{(^IscjgCls + UJrsclsCrs), (3) 
s 

Hqpc = ^(a;sfc4fcCsfc + i^Dqc'^gCoq) (4) 

kq 

are, respectively, the free Hamiltonians of the DQD, 
the electrodes coupled to the DQD, and the QPC with- 
out the tunneling term. In the DQD Hanhltonian, e is 
the energy detuning between the two single-dot levels 
and n is the intcrdot tunneling-coupling strength. Also, 



we define pseudospin operators = 0.20,2 — a\oi and 
Gx = a\ai + a\a2, with oi {02) being the annihilation 
operator for an electron staying at the left (right) dot. 
c-is (crs) is the annihilation operator for electrons in the 
source (drain) reservoir, i.e., the left (right) electrode of 
the DQD, while csk (coq) is the annihilation operator for 
electrons in the source (drain) reservoir of the QPC with 
momentum k (g). The tunneling-coupling Hamiltonian 
between the DQD and the two electrodes is 

Ht = Y^misclai + nrsTlcla2) + H.c], (5) 

s 

where the counting operator (Tj) decreases (in- 
creases) the number of electrons that have tunneled into 
the right electrode (via the barrier between the DQD and 
the right electrode).— They are introduced to keep track 
of the progress of the tunneling processes by successive 
electrons. Finally, 

Hdct = - Cct^)(4/cCD9 + C^qCsfc), (6) 

kq 

describes tunnelings in the QPC. It depends on the 
electron occupation of the DQD, owing to electrostatic 
couplings between the DQD and the QPC. We define 
T = To - (C2 + Ci)/2 and C = (C2 - Ci)/2, so that the 
transition amplitudes of the QPC, when an extra elec- 
tron staying at the left and right dots, equals T + ^ and 
T — C, respectively!^ 

To study the FCS of the transport through the DQD 
system, it is essential to know the probability P{n, t) of 
n electrons having been transported from the DQD to 
the right electrode during a period of time t. It is related 
to the cumulant generating function (CGF) G(x) defined 

e-G(x) =^p(„^i)eixn^ (7) 

n 

We will consider the time interval t much longer than 
the tunneling time of an electron through the DQD 
system, so that transient properties are insignificant. 
The derivative of the CGF with respect to the count- 
ing field X at X = yields the j-th cumulant, i.e. 
Cj = — (— i(9x)^G(Y)|x_>o, where x is a field conjugate to 
n (see, e.g., Ref. y)- These cumulants carry complete in- 
formation on the FCS of the DQD system. For instance, 
the average current and the shot noise can be expressed 
as / = eCi/t and S ~ 2e^C2/i. Thus, the Fano factor 
F, which is used to characterize the bunching and anti- 
bunching phenomena in the transport processes, is given 
hy F = S/2eI = C2/C1. The asymmetry (or skewness) 
of the distribution is given by the third-order cumulant 
C3. 

On the other hand, the probability distribution func- 
tion of the number of transported charges can also be 
expressed as 

P(n,t) = p^^o(t)+p;'^(t)+p:,(t), (8) 
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where («, j € {0,g,e}) denote the reduced density we can convert the master equation to 

matrix elements of the DQD at a given number n of elec- 
trons transported from the DQD to the right electrode 
in time t. Here 0, g, and e denote the eigenstates |0), 
and |e) of the DQD, which corresponds to no electron 

staying in the DQD, one electron in the ground state, „ 

and one electron in the excited state, respectively [see — = —Ai{x)Q, (10) 



horizontal solid lines in Fig. [IJb)] . These reduced den- 
sity matrix elements satisfy a master equation derived in 
Appendix A. Fourier transforming the reduced density 
matrix elements using 

ft,(x,i)-^pl}(0e^^", (9) with 



Mix) - -/JTi -7,, 




-a^Ffle^^f 

-Ire 
Ct'^^R + Ire 





2a/3Ffle^'^ 
-afiVn - 2777de 



7do 



A 



\Vr 4- 277^7^^ - 7^^ - 7j.^ / 



(11) 



where q = (poo(x, *), Pss(x, *), Pee(x, t), Re[peg(x, 0], 
Im[/9eg(X7 ^)])"^- In the matrix A^(x), ol = cos(6'/2), 
/3 = sin(6'/2), and 77 = cosS, with tan6' = 2i7/£; 
A = Ve^ -I- 4ri2^ and rL(fl) = '^'^9Li^B)^]^k{rk) ^^e rate 
of electron tunneling through the barrier between the 
DQD and the left (right) electrode. Here, gi (i = L, 
or R) denotes the density of states at the left or right 
electrode of the DQD, which is assumed to be constant 
over the relevant energy range. The QPC-induced exci- 
tation rate 7ex, relaxation rate 7rc, and dephasing rate 
7do are given by 



7ex = A [e (eVQPc - A) 4 
7rc = A [e (eVqpc + A) 4 
7dc = A[e (eFQPc) + e(- 



e (-e^QPc - A)] 
e (-e^QPc + A)] 
-eVQPc)] , 



(12) 



where A — 2TrgsgDC^ , with gs (go) being the density of 
states of the source (drain) electrode in the QPC. We 
also assume gs and gn to be constant over the relevant 
energies.^ Finally, Q{x) = (x + \x\)/2. 

The dynamical evolution of g(xy t) at long time t is 
governed hy g = -Amin(x)£'i where Amin(x) (> 0) is the 
minimal eigenvalue of A4{x)- Thus, 



tXn 



POoiX, t) + PggiX, t) + Peeix, t) 



Ce~ 



»(x)t 



(13) 



Since probability normalization implies P{n^ to) = Ij 
we have C = 1 and Amin as x — > 0. Then, from 
Eq. ([7]), the generating function G(x) can be obtained 
from^ G(x) = A^ini- 

With the transition rate of a single electron hopping 
from one reservoir of the QPC to the other, we can use 



the Landauer formula to obtain the transition probabil- 
ity^ and then have gsgo ~ -^QPC^/27rT^e^VQpc, where 
-^QPC and Vqpc arc the current and bias voltage of the 
QPC with densities of states gs and at the source and 
drain reservoirs. Following a recent experiment reported 
in Ref. HO, we take 7= 500 uA and Vqpc = 0.5 meV. In 
addition, we choose T = 0.5 and C, = 0.022, as in Ref. [23l. 
so that the QPC conductance changes by ^ 1% if the 
number of electrons in the DQD changes by one<^ 



III. CHARGE-DETECTOR-INDUCED 
BACKACTION 

In order to obtain compact analytical results for the 
FCS, we consider, for simplicity, the case where energy 
difference between two single-dot levels is zero (i.e., e = 
0). For instance, the charge current through the DQD is 
obtained as 



and shot noise is 



S 



"3 ' 



(14) 



(15) 



where 

s = TlTI + 4 (2rL + Tr) - 2r lTr (7^^ + 7re) , (16) 

X (7cx + 7ro) [-r^j - 41)2 + Tr (7ex + 7rc)] ■ (17) 
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Then, the Fano factor F = S/2eI also follows straight- 
forwardly. 

From Eqs. (14)-(17), it is clear that the charge current 
/, the shot noise 5", and hence also the Fano factor F 
depend on the excitation rate 7ex and the relaxation rate 
7rc induced by the QPC. These reveal the impacts of the 
backaction from the charge detector. More importantly, 
because of the nontrivial dependence of both 7ox and 7ro 
on the applied voltage across the QPC, the presence of 
the charge-detector-induced backaction could be experi- 
mentally checked more easily. This will be explained in 
detailed in the following. In addition, as a comparison, 
we also consider the no-backaction cases in which 7ox and 
7re in Eqs. ([T^ and (fTS)) equal zero. These no-backaction 
results are in accord with previous results given by other 
approaches j^^^ In our calculations, parameters like the 
interdot coupling strength and the tunneling rate Fl 
are taken from the experimental data whenever possi- 
ble3^ 

The charge current obtained from the cumulant Ci of 
the FCS is calculated both with and without backaction 
and the results are presented in Fig. [2] When the back- 
action from the charge detector is taken into account, we 
observe that the current through the DQD is significantly 
enhanced. In particular, when leVqpcl < A, a plateau 
with a constant current is observed. This plateau corre- 
sponds to a regime in which QPC-induced excitations is 
suppressed while there is still a constant relaxation rate 
contributed by the presence of the QPC, as can be in- 
terpreted from Eq. p2)) . Physically, a critical energy A 
exists for the QPC-induced excitation of an electron in 
the DQD and is hence required to change the current i ^^i^'^ 
Beyond the constant regime, that is, leVqpcl > A, it is 
clearly shown that the current increases with the magni- 
tude of the voltage applied across the QPC. 

Compared with the case of Tji = O.OSmeV [see Fig. 
[2Ua)], if the rate of electron tunneling to the right elec- 
trode is tuned to larger values, e.g., F/j = 0.15meV 
[Fig. IHb)] and 0.25meV [Fig. [^fc)], the current enhance- 
ment due to the backaction becomes more pronounced. 
This behavior may be tractable by directly examining 
Eq. Physically, it can be understood since an es- 

sential precondition for current flow in the unidirectional 
electron transport [see Fig. HJa)] is dominated by the in- 
teraction between the DQD and the right electrode. 

For the Fano factor F = S/2eI, our findings are given 
in Fig. [3l Without backaction, e.g., when the QPC is 
decoupled to the DQD, the Fano factor is always smaller 
than 1 (sec dotted lines in Fig. [3]). This means that anti- 
bunching of electrons corresponding to sitfe-Poissonian 
noise is expected (the physical reasons will be discussed 
below). If QPC-induced backaction is considered but 
with a condition [eVqpcl < A, we find a plateaux similar 
to that in the current. As mentioned above, the electron 
transport in this regime does not involve QPC-induced 
excitations. Beyond the threshold (i.e. leVqpcl > A), 
the Fano factor increases with |Vqpc|- For a sufficiently 
large bias, we can get _F = 1 [see Figs, ^h) and^c)]. 
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FIG. 2: (Color online). Current / versus the QPC bias energy 
eVqpc for tunneling rate r_R equal to: (a) 0.05 meV, (b) 0.15 
meV, (c) 0.25 meV. We use typical experimental parameters 
n = 0.1 meV and Tl = 0.05 meV from Refs. [H and [2l|. 



indicating that the electron transport is uncorrelated in 
time and is described by Poissonian statistics. Beyond 
this second threshold, F > 1. Thus, bunching of elec- 
trons in the transport through the DQD occurs, resulting 
in super-Poissonian noise. This QPC-induced change of 
shot noise from being sub- Poissonian (i^ < 1) to super- 
Poissonian {F > 1), and vice versa, becomes more pro- 
nounced with a smaller threshold voltage when the tun- 
neling rate F/^ is increased. 

Besides the charge-detector-induced backaction dis- 
cussed above, one might expect that dynamical chan- 
nel blockade could also be responsible for the super- 
Poissonian noise. In the following we rule out this pos- 
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FIG. 3: (Color online). Fano factor F versus QPC bias en- 
ergy eVqpc for tunneling rate Fr equal to: (a) O.OSmeV, (b) 
0.09 meV, (c) 0.15 meV. Parameters are the same as those in 
Fig. [3 



However, if the bias voltage applied across the QPC in- 
creases, transitions between the two eigenstate channels 
are enhanced [see Eq. ([T^ ] since 7ox = A|eVQPc — A| 
and 7rc = A leVqpc + A| when |eVQpc| > A . Then an 
electron can always switch to the faster tunneling chan- 
nel. For example, if \g) has a significantly higher tun- 
neling rate, the transport will mainly proceed through 
|0) ^ |e) ^ \g)^^ |0) and |0) ^ \g) ^ |0). ^ That 
means that two eigen-state channels are not equivalent 
due to different effective tunneling rates. Thus, super- 
Poissonian noise as shown in Fig. [3] is an expected result, 
fnterestingly, this backaction-induced changes of electron 
transport between bunching and anti-bunching behaviors 
should be realizable in the experiments with current tech- 
nologies. 

Note that the QPC-induced dephasing is not discussed 
above. This is because such dephasing with a rate 7de 
AleVqpcl) docs not induce any transitions between 
different states of the DQD. Finally, we emphasize that, 
besides using an approach different from the previous 
ones (see Appendix), our results are derived from the 
FCS. The formalism can be extended to the study of high 
order correlations. We can also go beyond zero-frequency 
counting statistics presented above and further study 
frequency-dependent current statistics in the presence of 
charge-dctector-induccd backaction. This will fully char- 
acterize the fundamental tunneling processes including 
transient behaviors and will be discussed elsewhere. 



IV. SUMMARY 

We have studied the unavoidable charge-detector- 
induced backaction on the FCS of a biased DQD. We 
find that this backaction has profound effects on the 
FCS, namely enhancing the current through the DQD 
and changing the shot noise from a sub-Poissonian regime 
to a super-Poissonian one. These backaction effects can 
be experimentally examined by using the current tech- 
nologies. Also, our results contribute to possible fine 
manipulation of quantum transport processes. When the 
measurement device is absent, we recover known results 
for both current and shot noise. 



sibility and emphasize that the behaviors of the shot 
noise observed above are due to the QPC-induced back- 
action. When the QPC is decoupled to the DQD, the 
rates of electron tunneling to the right electrode from 
the two eigen-state channels \g) and |e) in the DQD are 
given from the dynamical equation (see Appendix) as 

F^^ = P^Tr and F^^^ = a^Tu with a ^ ^, in gen- 
eral, which will lead to swper-Poissonian noise due to dy- 
namically blocked channelsi^"— However, as emphasized 
above, the assumption of a zero energy difference between 
two orbital levels (i.e., e = 0) leads to a = /3. This leads 
to sttfe-Poissonian noise. This is exactly what we have 
seen in Fig. |3] when the backaction is small or absent. 
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Appendix A: Quantum dynamics of the DQD 

Here, we derive a master equation to describe the quan- 
tum dynamics of the DQD. which is used to calculate the 
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FCS. In the eigenstate basis, the DQD HamHtonian can 
be written as 

HBQB = j{\e){e\~\g){g\), (Al) 

where A ~ \/e^ + 4rP is the energy sphtting of the two 
eigenstates of the DQD given by \g) = a\l) — I3\2), and 
|e) = (3\1) + a\2). In the interaction picture with the un- 
perturbed Hamihonian Hq = i?DQD + -ff leads + -ffQPCi the 
interaction Hamihonian Hi = iJx + -ffdct can be written 
as 

ifdct =X{t)Y (t) , (A2) 

S 

X (aaee*^*/2 - /3age-*^*/2*)e"""* + H.c] ,(A3) 

where 

3 

Tl=l 

kq 

and L/i =C|e)(g|,f/2 = Clff) C/3 - T - C cos 0^, = 
= A,W3 = 0, (t) = c^o^cske-'^^^"-^""^* . 
Applying the Born-Markov approximation and tracing 
over the degrees of freedom of the QPC, the quantum 
dynamics of the DQD system in the Schrodinger picture 
is governed by the master equation, 

p (t) = -^ [i/DQD ,p{t)]+ CdP (t) + CTp{t) , ( A6) 

with 

3 

^dpW- E {^^[P^]p{t) + ^^[P^,P3]p{t)} 

xA [e (eVqpc - uji) + 6 (-eVgpc - uji)] , 

CTp{t) = a^rLV[al]p {t) + l3''rLV[al]p (t) 

+a^TRV[a,T^p {t) + P^TRV[agT^p (t) 

+af3rL{ [al,p{t)ag] + [alp{t),ae] 

+ [al,p{t)a,] + [alp{t),ag] } 

-a/3rfl{ [aeTlp{t)alT] + [a^T^ p {t) , alT] 

+ [alT\p{t)alT] + [a.T^p(t),alT]}. 



Here, p (t) is the reduced density matrix of the DQD 
system. The superoperator I?, acting on any single or 
double operator, is defined as 

V[A]p^ApA^ -^A^Ap-^pAU, (A7) 
V[A,B]p=^ {ApB^ + BpA^ -B^Ap- pA^B) . (A8) 

From Eq. (|A6|) and the relations 

{n\TlpTr\n) - p^"-^\ (n|T,pTt|n) = p("+i),(A9) 
{n\TlTrp\n) = p^, (n|T,Ttp|n) = (AlO) 

where n is the number of electrons that have tunneled to 
the right electrode, we obtain the n-resolved equation of 
motion for each reduced density matrix element: 

-a/3rfl(p(7i)+p(r''), (All) 

+ {laprn + r,jMf+p';;J), (A12) 

pif = P'Tlp^I^ + j,,p<y> - {a'Tn + lre)pif 

+ {laPTR~Tj^,,){p(j^^+p<y>), (A13) 

pif = -lApif + aPTLP^,l^ + (ia/3r« + vie^^J 

-l{le.+lre){pif~P<^J'J). (Al4) 

In particular, we obtain the tunneling rates for two eigen- 
state channels F^^^ = a2Fi,F^'^^ = I3^Tl,T^jC = P^Tr, 
and F^^ = a^Ffl. 



^ Y. Blanter and M. Biittiker, Phys. Rep. 336, 1 (2000). ^ Y. V. Nazarov (ed.). Quantum Noise in Mesoscopic 



7 



Physics (NATO Science Series, Kluwer, Dordrecht, 2003). 
^ L. S. Levitov and G. B. Lesovik, JETP Lett. 58, 230 
(1993). 

* L. S. Levitov, H. W. Lee, and G. B. Lesovik, J. Math. 
Phys. 37, 4845 (1996). 

^ W. Lu, Z. Ji, L. Pfeiffer, K. W. West, and A. J. Rimberg, 

Nature (London) 423, 422 (2003). 
® T. Fujisawa, T. Hayashi, Y. Hirayama, H. D. Cheong, and 

Y. H. Jeong, Appl. Phys. Lett. 84, 2343 (2004). 

J. Bylander, T. Duty, and P. Delsing, Nature (London) 

434,' 361 (2005). 

* S. Gustavsson, R. Leturcq, B. Simovic, R. Schleser, T. Ihn, 
P. Studerus, K. Ensslin, D. C. Driscoll, and A. C. Gossard, 
Phys. Rev. Lett. 96, 076605 (2006). 

^ C. Fhndt, C. Fricke, F. Hohls, T. Novotny, K. Netocny, T. 
Brandes, and R. J. Hang, Proc. Natl Acad. Sci. USA 106, 
10119 (2009). 

^° J. Gabelh and B. Reulet, Phys. Rev. B 80, 161203(R) 

(2009) . 

C. Fricke, F. Hohls, N. Sethubalasubramanian, L. Fricke, 
and R. J. Hang C Fricke, AppL Phys. Lett. 96, 202103 

(2010) . 

T. Choi, T. Ihn, S. Schon, and K. Ensshn, Appl. Phys. 
Lett. 100, 072110 (2012). 
" N. Ubbelohde, C. Fricke, C. Fhndt, F. Hohls, and R. J. 
Hang, Nature Communications 3, 612 (2012). 

D. A. Ivanov and A. G. Abanov, arXiv:1203.6325v2. 

Y. V. Nazarov and M. Kindermann, Eur. Phys. J. B 35, 
413 (2003). 

^® E. V. Sukhorukov, A. N. Jordan, S. Gustavsson, R. 

Leturcq, T. Ihn, and K. Ensslin, Nat. Phys. 3, 243 (2007). 
1^ H. O. Li, M. Xiao, G. Cao, C. Zhou, R. N. Shang, T. Tu, 

G. C. Guo, H. W. Jiang, and GuoPing Guo, Appl. Phys. 



Lett. 100, 092112 (2012). 

W. G. van der Wiel, S. D. Franceschi, J. M. Elzerman, T. 
Fujisawa, S. Tarucha, and L. P. Kouwenhoven, Rev. Mod. 
Phys. 75, 1 (2003). 

G. Kiefiich, P. Samuelsson, A. Wacker, and E. Scholl, Phys. 

Rev. B 73, 033312 (2006). 
^° P. Barthold, F. Hohls, N. Maire, K. Pierz, and R. J. Haug, 

Phys. Rev. L 96, 246804 (2006). 
2^ G. Kiefiich, E. Scholl, T. Brandes, F. Hohls, and R. J. 

Haug, Phys. Rev. Lett. 99, 206602 (2007). 

S. H. Ouyang, C. H. Lam, and J. Q. You, Phys. Rev. B 

81, 075301 (2010). 

S. Gustavsson et al, Phys. Rev. Lett. 99, 206804 (2007). 
S. A. Gurvitz and Y. S. Prager, Phys. Rev. B 53, 15932 
(1996). 

Y. V. Nazarov and J. J. R. Struben, Phys. Rev. B 53, 
15466 (1996). 

C. B. Doiron, B. Trauzettel, and C. Bruder, Phys. Rev. B 
76, 195312 (2007). 

Z. Z. Li, S. H. Ouyang, C. H. Lam, and J. Q. You, Phys. 
Rev. B 85, 235420 (2012). 

T. M. Stace and S. D. Barrett, Phys. Rev. Lett. 92, 136802 
(2004). 

S. A. Gurvitz, Phys. Rev. B 56, 15215 (1997). 
^° L. M. K. Vandersypen et ai, Appl. Phys. Lett. 85, 4394 
(2004). 

D. A. Bagrets and Yu. V. Nazarov, Phys. Rev. B 67, 
085316 (2003). 

A. Cottet, W. Belzig, and C. Bruder, Phys. Rev. Lett. 92, 
206801 (2004). 

A. Cottet and W. Belzig, Europhys. Lett. 66, 405 (2004). 
W. Belzig, Phys. Rev. B 71, 161301 (2005). 



